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In this article we study the dynamics of fermions in a fermion-scalar plasma. We begin by 
obtaining the effective in-medium Dirac equation in real time which is fully renormalized and causal 
and leads to the initial value problem. For a heavy scalar we find the novel result that the decay 
of the scalar into fermion pairs in the medium leads to damping of the fermionic excitations and 
, their in-medium propagation as quasiparticles. That is, the fermions acquire a width due to the 

■ decay of the heavier scalar in the medium. We find the damping rate to lowest order in the Yukawa 

^\ ' coupling for arbitrary values of scalar and fermion masses, temperature and fermion momentum. 

t-H , An all-order expression for the damping rate in terms of the exact quasiparticle wave functions is 

(— i ■ established. A kinetic Boltzmann approach to the relaxation of the fermionic distribution function 

' confirms the damping of fermionic excitations as a consequence of the induced decay of heavy 

scalars in the medium. A linearization of the Boltzmann equation near equilibrium clearly displays 
' the relationship between the damping rate of fermionic mean fields and the fermion interaction rate 

H , to lowest order in the Yukawa coupling directly in real time. 

(SJ ■ 12.15Ji,12.38.Mh 
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1 I. INTRODUCTION 

o : 

qq , The propagation of quarks and leptons in a medium of high temperature and/or density is of fundamental im- 
• portance in a wide variety of physically relevant situations. In stellar astrophysics, electrons and neutrinos play a 
major role in the evolution of dense stars such as white dwarfs, neutron stars and supernovae jjj. In ultrarelativistic 
heavy ion collisions and the possibility of formation of a quark-gluon plasma, electrons (and muons) play a very 
important role as clean probes of the early hot stage of the plasma Q . Furthermore, medium effects can enhance 
neutrino oscillations as envisaged in the Mikheyev-Smirnov-Wolfenstein (MSW) effect || and dramatically modify 
the neutrino electromagnetic couplings Q . The propagation of quarks during the non-equilibrium stages of the elec- 
troweak phase transition is conjectured to be an essential ingredient for baryogenesis at the electroweak scale both in 
non-supersymmetric and supersymmetric extensions of the standard model |5| . 

In-medium propagation is dramatically different from that in vacuum. The medium modifies the dispersion relation 
of the excitations and introduces a width to the propagating excitation |3|-^j that results in damping of the amplitude 
of the propagating mode. In this article we focus on several aspects of propagation of fermionic excitations in a 
fermion-scalar plasma: 

• We begin by deriving the effective and fully renormalized Dirac equation in real time. This is achieved by 
relating the expectation value of a fermionic field induced by an external fermionic source via linear response 
to an initial value problem for the expectation value. This initial value problem is in terms of the effective real 
time Dirac equation in the medium that is i) renormalized, ii) retarded and causal. 

The necessity for a consistent Dirac equation in a medium has been recognized in the literature within the 
context of neutrino oscillations in the medium |Io|-|l3[|. In particular, in references [^o|JTl| l a proposal for a field- 
theoretical treatment of neutrino oscillations in the medium starting from the Dirac equation was presented. 

In this article we introduce the fully renormalized, in-medium effective Dirac equation in real time that allows a 
more transparent study of damping and oscillations in time. There are definite advantages in such formulation 
since the real time evolution is obtained at once from an initial value problem and allows a straightforward 
identification of the damping rate. 

• We apply the effective Dirac equation in a medium to study the real-time evolution of fermionic excitations 
in a fermion-scalar plasma. Whereas the propagation of quarks and leptons in a QED or QCD plasma has 
been studied thoroughly (see ref. MM for details) a similar study for a scalar plasma has not been carried 
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out to the same level of detail. Recently some attention has been given to understanding the thermalization 
time scales of bosonic and fermionic excitations in a plasma of gauge |14[| and scalar bosons |15| , furthermore 



fermion thermalization is an important ingredient in models of baryogenesis mediated by scalars |16|. Most 
of the studies of fermion thermalization focus on the mechanism of fermion scattering off the gauge quanta in 
the heat bath and/or Landau damping in the hard thermal loop (HTL) resummation program JL7J . Although 
the scalar contribution to the fermionic self-energy to one loop has been obtained a long time ago 0, scant 
attention has been paid to a more detailed understanding of the contribution from the scalar degrees of freedom 
to the fermion relaxation and thermalization. As mentioned above this issue becomes of pressing importance in 
models of baryogenesis and more so in models in which the scalars carry baryon number Jl6| . 

Whereas the contribution from scalars to the fermionic thermalization time scale (damping rate) has been studied 
for massless chiral fermions JTojJT^ ] , in this article we offer a detailed and general study of fermion relaxation 
and thermalization through the interactions with the scalars in the plasma in real time and for arbitrary values 
of the scalar and fermion masses, temperature and fermion momentum. More importantly, we focus on a novel 
mechanism of damping of fermionic excitations that occurs whenever the effective mass of the scalar particle 
allows its kinematic decay into fermion pairs. This phenomenon only occurs in a medium and is interpreted as 
an induced decay of the scalar in the medium. It is a process different from collisions with particles in the bath 
and Landau damping which are the most common processes that lead to relaxation and thermalization. 

This process results in new thermal cuts in the fermionic self-energy and for heavy scalars, this cut results in a 
quasiparticle pole structure for the fermion and provides a width for the fermionic quasiparticle. The remarkable 
and perhaps non- intuitive aspect of this process is that the decay of the scalar results in damping of the fermionic 
excitations and their propagation as quasiparticle resonances. Our real time analysis reveals that amplitude of 
the /c-mode of the expectation value of the fermion field decreases as e~ Tk * while it oscillates with frequency 
Ld p (k) which is determined by the position of the resonance. 

The effective real-time Dirac equation in the medium allows a direct interpretation of the damping of the quasi- 
particle fermionic excitation and leads to a clear definition of the damping rate. By analyzing the quasiparticle 
wave functions we obtain an all-order expression for the damping rate Tk and confirm and generalize recent 
results for the massless chiral case | ]To| , ^8| . 

• In order to provide a complementary understanding of the process of induced decay of the heavy scalars in the 
medium and the resulting fermion damping, we study the kinetics of relaxation of the fermionic distribution 
function via a Boltzmann equation to lowest order. Linearizing the Boltzmann equation near the equilibrium 
distribution, we obtain the relation between the thermalization rate for the distribution function in the relaxation 
time approximation (linearized near equilibrium) and the damping rate for the amplitude of the fermionic mean 
fields to lowest order in the Yukawa coupling. This analysis provides a real-time confirmation of the oft quoted 
relation between the interaction rate (obtained from the Boltzmann kinetic equation in the relaxation time 
approximation) and the damping rate for the mean field More importantly, this analysis reveals directly, 

via a kinetic approach in real time how the process of induced decay of a heavy scalar in the medium results in 
damping and thermalization of the fermionic excitations. A study of the relation between the interaction rate 
and the damping rate has been presented recently for gauge theories within the context of the imaginary time 
formulation Jl{J. Our results provide a real-time confirmation of those of reference jl9| for the scalar case. 

The article is organized as follows: in section II we obtain the effective in-medium Dirac equation in real time 
starting from the linear response to an external Grassmann-valued source that induces a mean field. We obtain the 
fully renormalized Dirac equation with the real time self-energy to one loop order by turning the linear response 
problem into an initial value problem for the mean field. The renormalization aspects are addressed in detail in this 
section. In section III we study in detail the structure of the renormalized self-energy and establish the presence of new 
cuts of thermal origin. We then note that for heavy scalars such that their decay into fermion pairs is kinematically 
allowed, the fermionic pole becomes embedded in this thermal cut resulting in a quasiparticle (resonance) structure, 
which is analyzed in detail. The decay rate is evaluated in the narrow width approximation (justified for small Yukawa 
couplings) for arbitrary values of the scalar and fermion masses, temperature and fermion momentum. 

In section IV we present a real time analysis of the evolution of the mean- fields. In this section we clarify the 
difference between complex poles and resonances (often misunderstood) . This analysis reveals clearly that the induced 
decay of the scalar results in an exponential damping of the amplitude of the mean field and yields to a clear 
identification of the damping rate bypassing the conflicting definitions of the damping rate offered in the literature. 
An analysis of the structure of the self-energy and an interpretation of the exact quasiparticle spinor wave functions 
allows us to provide an all-order expression for the damping rate of the fermionic mean fields. In section V we present 
an analysis of the evolution of the distribution functions in real time by obtaining a Boltzmann kinetic equation 
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for the spin-averaged fermionic distribution function. In the linearized approximation near equilibrium (relaxation 
time approximation), we clarify to lowest order the relation between the damping rate of the quasiparticle fermionic 
excitations and the thermalization rate (linearized) of the distribution function, directly in real time. In the conclusions 
we summarize our results and suggest new avenues including the contribution from gauge fields. In this section we 
also assess the importance of the scalar contribution in theories with gauge and scalar fields. 

II. EFFECTIVE DIRAC EQUATION IN THE MEDIUM 

As mentioned in the introduction, whereas the damping of collective and quasiparticle excitations via the interactions 
with gauge bosons in the medium has been the focus of most attention, understanding of the influence of scalars has 
not been pursued so vigorously. 

Although we are ultimately interested in studying the damping of fermionic excitations in a plasma with scalars 
and gauge fields within the realm of electroweak baryogenesis in either the Standard Model or generalizations thereof, 
we will begin by considering only the coupling of a massive Dirac fermion to a scalar via a simple Yukawa interaction. 

The model dependent generalizations of the Yukawa couplings to particular cases will differ quantitatively in the 
details of the group structure but the qualitative features of the effective Dirac equation in the medium as well as the 
kinematics of the thermal cuts that lead to damping of the fermionic excitations will be rather general. 

We consider a Dirac fermion with the bare mass M coupled to a scalar with the bare mass uiq via a Yukawa 
coupling. The bare fermion mass could be the result of spontaneous symmetry breaking in the scalar sector, but for 
the purposes of our studies we need not specify its origin. 

The Lagrangian density is given by 

C = - M Q )<f + \d^d^ - ^ml4> 2 - CM - y *<l>* + ^ + *r, + j<f> , (2.1) 

where yo is the bare Yukawa coupling. The self-interaction of the scalar field accounted for by the term Ci[<fi] need not 
be specified to lowest order. The r\ and j are the respective external fermionic and scalar sources that are introduced 
in order to provide an initial value problem for the effective Dirac equation. We now write the bare fields and sources 
*, 4>, r) and j in terms of the renormalized quantities (hereafter referred to with a subscript r) by introducing the 
renormalization constants and counterterms: 

* = Z)l 2 * r , = z\ /2 <t> r ,rj = Z- 1/2 Vr , j = Z~ 1/2 Jr , 

y = y zl /2 Z^/Z y , m 2 = (S m + m 2 ) /Z^ , M = (S M + M) /Z* . (2.2) 
With the above definitions, the C can be expressed as: 

C = 9 r (i@- M)* r + ^d^r d^<f) r - ^m 2 (f)l - C r j[<j> r } - y* r cj) r * r + fj r * r + * r r] r + j r <t> r 

+-<^d M r <3^0 r - -<5 m <^ + ^ , (2.3) 

where m and M are the renormalized masses, and y is the renormalized Yukawa coupling. The terms with the 
coefficients 

S,p — Z^p — 1 , 5$ — Ztf, — 1 , 

S M = MqZ^ - M , S m = mlZ^ - m 2 , (2.4) 
$y = z y ~ 1 

and SUj are the counterterms to be determined consistently in the perturbative expansion by choosing a renormal- 
ization prescription. As it will become clear below this is the most natural manner for obtaining a fully renormalized 
Dirac equation in a perturbative expansion. 

The dynamics of expectation values and correlation functions of the quantum field is obtained by implementing 
the Schwinger-Keldysh closed-time-path formulation of non-equilibrium quantum field theory p0|-p3|. The main 
ingredient in this formulation is the real time evolution of an initially prepared density matrix and its path integral 
representation. It requires a path integral defined along a closed time path contour. This formulation has been 
described elsewhere within many different contexts and we refer the reader to the literature for details [po|-p3| . 
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Our goal is to understand the non-equilibrium relaxational dynamics of the inhomogeneous fermionic mean fields 

l/>(2,t) = (*r(f,t)) 

from their initial states in the presence of the fermion-scalar medium. 

This statement requires clarification. In states of definite fermion number (either zero or finite temperature) the 
expectation value of the fermion field must necessarily vanish. However, in order to understand the non-equilibrium 
dynamics, we prepare the system by coupling an external Grassman source to the fermionic field in the Hamiltonian. 
This source creates a coherent state of fermions which is a superposition of states with different fermion number. As 
will be discussed in detail below, this source term is switched-on adiabatically from t = — oo in the full interacting 
theory up to a time to and switched-off at t — 4 . The resulting fermionic coherent state at t — 4 does not have a 
definite fermion number and the expectation value of the field operator in this state is non- vanishing. After the source 
is switched-off at time to the expectation value evolves in time with the full interacting Hamiltonian in absence of 
sources. This is the standard approach to linear response, in this approach the source term couples to a given mode 
of wavevector k of the fermionic fields, thus displacing this degree of freedom off-equilibrium. The other modes are 
assumed to be remain in thermal equilibrium. 

Since we are interested in real time correlation functions and the initial density matrix is assumed to be that of 
thermal equilibrium at initial temperature T — with respect to the free (quadratic) Lagrangian, only the real 
time branches, forward and backwards are required. The contribution from the imaginary time branch corresponding 
to the thermal component of the density matrix cancels in the connected non-equilibrium expectation values [pl|-p3"[ . 
The effective non-equilibrium Lagrangian density that enters in the contour path integral is therefore given by 

£„o„- oq = c [*+, §+, 0+] - c [*-, *-,#:] . (2.5) 

Fields with (+) and (— ) superscripts are defined respectively on the forward (+) and backwards (— ) time contours 
and are to be treated independently The external sources are the same for both branches. 

The essential ingredients for perturbative calculations are the following real time Green's functions (23) : 

• Scalar Propagators 



G t 


+(t,t')- 


= G|(4,4')6(4- 


4') + G|(4,4')6(4'- 
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(1 + n k )e 


-iw fe (t— t') _|_ nke ^k(t- 


-t') 




(4, t') = 




3 x e~ zts ( 


4>r(6, t')<j) r {x, t)) 








i 

2u k 




+ (1 + n k )e tuJk(t - 


-f) 



u> k - \jp + m i , n k - -^J— - . (2.6) 



• Fermionic Propagators (Zero chemical potential) 

S++ (4, t') = 5| (4, 4') 9(4 - 4') + Sf(t, t')<d(t' - t) , 
^ "(*>0 = ^(t,t')Q(t' -t) + S5(4,4')6(4 - t') , 
S+-(t,t') = -55(4,4'), 
SZ+(t,t') = -SZ(t,t'), 

51(4,4') = -i f d 3 x e- its (* r (x,4)* r (0,4')) 
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(Jf + M)(1 - n fc )e-* a *( t - t,) + 7o(# ~ M^e^*-*') 



I 

2Q k 



(ff + M)n k e- iQ ^-^ + 7 o(^ - M) 7o (l - r^y^*-*') 
~ 1 



"A- 



-|- 1 



(2.7) 



The perturbative evaluation of correlation functions proceeds as usual, but now the Feynman rules involves two 
types of vertices with opposite signs and the four different non-equilibrium propagators. The symmetry factors are 
the usual ones. 

These free propagators ( [2.6| ) and ( |2.7D are thermal since the initial state is assumed to be in thermal equilibrium. 



A. In-medium Dirac equation from linear response 

Consider the fermionic mean field obtained as the linear response to the externally applied (Grassmann-valued) 
source 7] r : 

/oo 
dt'd 3 x'S Ict {x - x', t - t 1 ) r, r (x', t') , (2.8) 
-oo 

with the exact retarded Green's function 

S Ict (x-x',t~t') = [S > (x-x',t-t')-S < {x-x',t-t')] &{t-t') 

= -i{{y r (x,t),y r (x , ,t , )})Q(t-t') , (2.9) 

where the expectation values are in the full interacting theory but with vanishing sources. An initial value problem 
is obtained by considering that the external fermionic sources are adiabatically switched on in time from t — > — oo 
thereby inducing an expectation value of the fermionic fields, and switching-off the source term at some time to. Then 
for t > to this expectation value or mean field will evolve in the absence of a source and will relax because of the 
interactions. The evolution for t > to is an initial value problem, since the source term was used to prepare an initial 
state and switched off to let this state evolve in time. This initial value problem can therefore be formulated by 
choosing the source term to be of the form 

ri r (x, t) = ? ]r (x) e et Q(t - t) r, e -> 0+ . (2.10) 

In what follows we choose to — for convenience. The adiabatic switching-on of the source induces an expectation 
value that is dressed adiabatically by the interaction. The retarded and the equilibrium nature of S Ie t(x — x' , t — t') 



(which depends on the time difference) and the form of the source (2.10) guarantee that 

i/j(x,t = 0) = ijj (x) , 
ip(x,t < 0) = , 

where ij)a{x) is determined by rj r (x) (or vice versa, the initial conditions for the condensates ipo{%) can be used to 



find rj r {x)). This can be seen by taking the time derivative of ip(x,t) in eq. (2.8), using the fact that the retarded 
propagator depends on the time difference, integrating by parts, and using the form of the external source and the 
retarded nature of the propagator. 

In order to relate this linear response problem to the initial value problem for the dynamical equation of the mean 
field, let us consider the (integro) differential operator 0^,t) which is the inverse of S le t(x — x', t — t') so that 

0( S ,tMx,t) = -r)r(x,t) , ip(x,t = 0) = i) Q {x) , ^(2,t<0) = 0, (2.11) 



where the source is given by eq. (2.10) 



It is at this stage where the non-equilibrium formulation provides the most powerful framework. The real-time 
equations of motion for the mean fields can be obtained via the tadpole method |2l|-|23|], which automatically leads 
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to a retarded and causal initial value problem for the expectation value of the field. The implementation of this 
method is as follows. Let us introduce the inhomogeneous mean fields ip(x,t) = (&(x, t)) and ip{x,t) = fi?(x, t)}. The 
dynamics of these fermionic mean fields in the plasma can be analyzed by treating ip(x,t) and 4>{x,t) as background 
fields, i.e., the expectation values of the corresponding fields in the non-equilibrium density matrix, by expanding 
the non-equilibrium Lagrangian density about these mean fields. Therefore, we write the full quantum fields as the 
c-number expectation values (mean fields) and quantum fluctuations about them: 



#± (x, t) = xP(x, t) + x ± (x, t) , ^± (x , t) = i){x, t) + x± (x, t) 
i/>(2, t) = (tf±(£, t)) , $(2, t) = (*±(z, t)) , 



(2.12) 



where the expectation values of the field operators are taken in the time-evolved density matrix. 

The equations of motion for the mean fields can be obtained to any order in the perturbative expansion by imposing 
the requirement that the expectation value of the quantum fluctuations in the time ev olved density matrix vanishes 
identically. This is referred to as the tadpole equation |2l|-f23| which follows from eqs. fl2.12 ): 



(x ± ) = o, (x ± ) = o. 



(2.13) 



The procedure consists in treating the linear terms in x as well as the non-linearities in perturbation theory and 
keeping only the connected irreducible diagrams, as in usual perturbation theory. The equations obtained via this 
procedure are the equations of motion obtained by variations of the non-equilibrium effective action |2^-|2^] . 

By applying the tadpole method and taking spatial Fourier transform, we find the effective real time Dirac equation 
for the fc-mode of the expectation value of the fermion: 



where — t') is the fermion self-energy and 



>M 



Mt) + / dt' ~ il>M) = - V r(t) , 



Mt)= / d s xe- %k - x i>(Z,t) . 



Using the non-equilibrium Green's functions (|2.6|) and (2.7), we find to one loop order, that Sr(i — f) is given by 



where 



Eg(i-0=Ho4 U; (t-i')+7-fe Y-fit-n + Vfit-t') , 



4 0) (*-o 



4 2) (*-<o 



V 



dh. 



(2n) 3 2uk+qUi q 



y 



COs[(uJk+q + LOq){t - <')](1 + 7l k+q - H q ) + COs[(uJ k+q - LU q )(t - t')](n k+q + fig) 

d 3 q 1 k ■ q 



V 



(27r) 3 2bJk+q&q k 2 

sm[(uj k+q + u> q )(t - t')](l + n k+q - n q ) - sm[{u k+q - Q q )(t - t')](n k+q + n q ) 
d 3 q M 

(27r) 3 2(Jk+q&q 

sm.[(u k +q + u) q ){t - t')](l + n k+q - n q ) - sin[(w fc+9 - Q q )(t - t')](n k+q + n q ) 



(2.14) 



with uj k + q = \J (k + q) 2 + m 2 and n k+q = (e^ UJk+q ~ 1) 1 being, respectively, the energy and the distribution function 

for scalars of momentum k + q. 

As mentioned before, the source is taken to be switched on adiabatically from t = — oo and switched off at t = to 
provide initial conditions 



= °) = M°) > < o) = o . 



(2.15) 



Defining a^(t — t') as 



G 



—<r S (t = - 

and imposing that ?7 r g(t > 0) = 0, we obtain the equation of motion for t > 



(2.16) 



i"(o-Qj_ - 1 ■ k - M j +8^ I *7ot^ - 7 • fc ) + ^(0) - cS M 



^(t) - / ag(t - t') ^ k (t') = 
Jo 



(2.17) 



The equation of motion (2.17) can be solved by Laplace transform as befits an initial value problem. The Laplace 
transformed equation of motion is given by 



ijos -j-k-M + S^ (ijos - 7 • - S M + erg (0) - sa%(s) ^g(s) 
= [ijo + iS^jo - erg(s)] ^g(O) , 

where i/>%(s) and er^(s) are the Laplace transforms of ^>n(i) and erg(t) respectively: 



(2.18) 



dt e~ st aAt) . 



B. Renormalization 



Before proceeding with the solution of the above equation, we address the issue of the renormalization by analyzing 
the ultraviolet divergences of the kernels. As usual the ultraviolet divergences are those of zero temperature field 
theory, since the finite temperature distribution functions are exponentially suppressed at large momenta. Therefore 
the ultraviol et di vergences are obtained by setting to zero the bosonic and fermionic occupation numbers. 

With eq. (12T61) , a % {t-t') can be written as 



a At -t') = l7 o4 0) (t-t')+i-k a^> (t-t')+ of (t - 



(i), 



J2), 



(2.19) 



A straightforward calculation leads to 

erf(0) = 0, a«(0) = -!§>(£) 



where crg(s), a^(s) and a¥'(s) are the Laplace transform of a y y'(t), ^'(t) and er^ J (t) respectively, A is an 
ultraviolet momentum cutoff and K, is an arbitrary renormalization scale. Therefore the counterterms 5^ and 5m are 
chosen to be given by 



finite 



f finite , crl 2) (0) = In ( A ) + finite , 
erl (s) = finite , o^ 2 ' 1 (s) — finite , 



.(2) 



16^ ln U 



finite , 5 



M 



V 2 M /A 



finite 



and the respective kernels are rendered finite, i.e. 



ct^(s) - iS^o = a r g(s) = finite , cr^(0) - ^7 • fc - cTat = er r ^(0) = finite 



(2.20) 



(2.21) 



The finite parts of the counterterms in eq. (2.2C) are fixed by prescribing a renormalization scheme. There are 
two important choices of counterterms: i) determining the counterterms from an on-shell condition, including finite 
temperature effects, and ii) determining the counterterms from a zero temperature on-shell condition. Obviously 
these choices only differ by finite quantities, however the second choice allows us to separate the dressing effects of the 
medium from those in the vacuum. For example by choosing to renormalize with the zero temperature counterterms 
on-shell, the pole in the particle propagator will have unit residue at zero temperature; however in the medium, the 
residue at the finite temperature poles (or the position of the resonances) are finite, smaller than one and determined 
solely by the properties of the medium. Thus the formulation of the initial value problem as presented here yields an 
unambiguous separation of the vacuum and in-medium renormalization effects. 

Hence we obtain the renormalized effective Dirac equation in the medium and its initial value problem 
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ilas - 7 ■ k - M + ^ rk {s) ^(s) 



«70 - <T r Hs 



Mo) , 



(2.22) 



with it „g(0) and cr g(s) the fully renormalized kernels (see eq. (|2.2l| )), and E r g(s) the Laplace transform of the 
renormalized fermion self-energy which can be written in its most general form as follows 



E rJf 00 - i 70 s ef (s) + 7 • k ef (s) + M ef (s) 
The solution to eq. ( 2.22 ) is 



l r 

s 



1 



-fc + M-I! rijE (0))]^(0) 



(2.23) 



(2.24) 



where 



S(s, k) = ij Q s - 7 • k - M + S n g(s) 

l7oS (i + 4 0) oo) - 7 • fca - 4 1} ( s )) + M (! - 4 2) ( s )) 



det 



(2.25) 



-s 2 (l + ef(s)) 2 - k 2 (l - ef(s)) 2 - M 2 (l - ef(s)) 2 
is the fermion propagator in terms of the Laplace variable s. The square of the denominator of ( |2.25| ) is recognized as 

ij s - 7 • k - M + S r 



The real-time evolution of ipj"(t) is obtained by performing the inverse Laplace transform along the Bromwich contour 
in the complex s-plane parallel to the imaginary axis and to the right of all singularities of ipj:(s). Therefore to obtain 
the real time evolution we must first understand the singularities of the Laplace transform in the complex s-plane. 



III. STRUCTURE OF THE SELF-ENERGY AND DAMPING PROCESSES 

To one loop order, the Laplace transform of the components if 1 (s) of the fermion self-energy Eg(s) (see eq. ( [2.23 )) 



fwi , j f 4 0) ( fc o) 

ef{s) = / dk Q - T —,\ k pf(k 



can be written as dispersion integrals in terms of spectral densities pi (ko) 



5a_ 

M 



with the one-loop spectral densities given by the expressions 

O i0) (kn)-V 2 [ d3<1 ^ X 

[<5(fc - w fc+9 - Q q )(l + n k+q - n q ) + S(k - uj k+q + Q q )(n k+q + n q )} 
(1),, , 2 f d 3 q 1 k-q 



* 7 (27r)- ;i 2w fe+9 o;g AH 

[<5(fc - u k+q - Q q )(l + n k+q - fig) -6(k - u k+q + ui q )(n k+q + n q )] 



4 a) (*o)=» a 



d 3 < 



* v u/ a 7 (2tt) 3 2cu k+q to q 

[S(k - uj k+q - LJ q )(l + n k+q - fig) - S(ko - to k+q + u)g)(n k +g + n q )] . (3.2) 

The analytic continuation of the self-energy and its components e^ (s) in the complex s-plane are given by 

^r,fc( S = ~ iu} ± 0+ ) = S fi^M ± i %I,k( U ) > 

sf (s = -iw ± 0+) = e «.( W ) ± * e W M , (3.3) 
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where the real parts are even functions of to and given by 



■2>> 

4m 



dfc 7> 



fc 2 _ W S 



I k pf(k ) 



M 



(3.4) 



and the imaginary parts are odd functions of to and given by 

4>) = ^M[4 0) (I-D + 4 0) H-I 



4 2 lM = ^sgnM [pf{\uj\)-pf{-\uj\) 



(3.5) 



The denominator of the analytically continued fermion propagator eq. (2.25) can be written in the compact form 

ui 2 - u>l + n(oj, k) 



with u>k = Vk 2 + M 2 and 



n( w , fc) = 2 oj 2 ef ( W ) + k 2 ef ( W ) + m 2 4 2) m 



+ c 2 [ef(c)] 2 - fc 2 [e«(c)] 2 - M 2 [ef (c)] 2 . (3.6) 
We recognize that the lowest order term of this effective self-energy can be written in the familiar form 



II(u/, k) = 2 k> 2 4°V) + fc a eg + M a ej^(w) I = ^Tr[(# + M)%(s = -iw + 0+)] 



.2 (i), 



r2„(2) 



(3.7) 



but certainly not the higher order terms. 

T he i maginary part of II on-shell will be identified with the damping rate (see belo w). The expression given by 
eq. (3.7) leads to the familiar form of the damping rate, but we point out that eq. (3.7) is a lowest order result. The 
full imaginary part must be obtained from the full function II(cj, k) and the generalization to all orders will be given 
in a later section below. 

For fixed M and m, the (S-function constraints in the spectral densities p2 (to) can only be satisfied for certain ranges 
of lu. Since the imaginary parts are odd functions of to we only consider the case of positive u>. The 6(\uj\ — LOk+q — & q ) 
has support only for \u>\ > \/k 2 + (m + M) 2 and corresponds to the normal two-particle cuts that are present at zero 
temperature corresponding to the process iji — ► 4> + ip- These cuts (both for positive and negative uj) do not give a 
contribution to the imaginary part on-shell for the fermion. 

The terms proportional to rik+q + n q give the following contribution to the lowest order effective self-energy: 



TLi(u),k) = 7ry 2 sgn(o;) 



d 3 q n k+q + n q 

(2tt) 3 2u) q LOk+q 



oj\uj q - k- q - M ) S (\uj\ - oj k+q + u> q ) 



+ ( \uj\uj q + k ■ q + M ) S (\uj\ + 0J k+q - iv q ) 



(3.8) 



The first delta function determines a cut in the region < |oj| < yk 2 + (m — M) 2 and originates in the physical 
process <j) — * ip + ip whereas the second delta function determines a cut in the region < \u>\ < k and originates in 
the process <fi + rp — > ip. Whereas the first cut originates in the process of decay of the scalar into fermion pairs, 
the second cut for ( lu 2 < fc 2 ) is associated with Landau damping. Both delta functions restrict the range of the 
integration variable q (see below). For m > 2M the scalar can decay on-shell into a fermion pair, and in this case the 
fermion pole is embedded in the cut < \u>\ < \Jk 2 + (m — M) 2 becoming a quasiparticle pole and only the first cut 
contributes to the quasiparticle width. This is a remarkable result, the fermions acquire a width through the induced 
decay of the scalar in the medium. This process only occurs in the medium (obviously vanishing at T = 0) and its 
origin is very different from cither collisional broadening or Landau damping. A complementary interpretation of the 
origin of this process as a medium induced decay of the scalars into fermions and the resulting quasiparticle width 
for the fermion excitation will be highlighted in section V within the kinetic approach to relaxation. 







The width is obtained to lowest order from Hi(uik,k), and the on-shell delta function is recognized as the energy 
conservation condition for the decay <f> — > %j> + tp of the heavy scalar on-shell. To lowest order we find the following 
expression for 11/ (uj k , k) for arbitrary scalar and fcrmion masses with m > 2M and arbitrary fermion momentum and 
temperature: 



ITj(wfc, k) = Try" 



d 3 q uj k Lu q — k-q — M 2 



(27r) 3 2u q CJ k+q 



(nk+q + n q ) S(uj k + Ul q ~ UJ k+q ) , 



y 2 m 2 T ( AM- 
167rfc 



In 



1 + e 



-fx 



(3.9) 



where q\ and q% are given by 



q 2 



m 

2M 2 

2 

777 

2M 2 



k 1 



k 1 




with \q\ S (q*, q%) being the support of 8{oj k — uik+q 



(3.10) 



IV. REAL TIME EVOLUTION 

The real time evolution is obtained by performing the inver se Laplace transform as explained in section II. This 
requires analyzing the singularities of ^r(s) given by eq. ( 2.25| ) in the complex s-plane. It is straightforward to see 
that the putative pole at s = has vanishing residue, therefore the singularities are those arising from the inverse 
fermion propagator S(s,k). If the fermionic pole is away from the multiparticle cuts, the singularities are: i) the 



isolated fermion poles at s 



with Lo p the position of the isolated (complex) poles (corresponding to stable 



fermionic excitations), and ii) the multiparticle cuts along the imaginary axis s = — iui for < \u>\ < \Jk 2 + (m — M) 2 
and M > y/k 2 + (m + M) 2 . 

The Laplace transform is performed by deforming the contour, circling the isolated poles and wrapping around the 
cuts. When the scalar particle can decay into fermion pairs, i.e., m > 2M, the fermion pole is embedded in the lower 
cut and we must find out if it becomes a complex pole in the physical sheet (the domain of integration) or moves off 
the physical sheet. 

For to < 2M the fermion pole at ui = uj p is real and the fermion mean field oscillates for late times with constant 
amplitude and frequency oj p . 



A. m > 2M: complex poles or resonances? 



When m > 2M the fermion pole is embedded in the cut < \u>\ < ^Jk 2 + (m — M) 2 and the pole becomes 



complex. The position of the complex poles are obtained from the zeros of u> 



II(u;, k) in the analytically 



continued fermion propagator for lj = 0J v (k) — iY k with lo p being the real part of the complex pole. For T k <C ld p 
(narrow width approximation) and with the expressions for the discontinuities in the physical sheet given by eq. ( \i.<\ ) , 
the equation that determines the position of the complex pole is given by the solution of the following equation 



(w p - i T k ) 2 - lo\ + n fl (w p , k) - i sgn(r fc ) Ui(lu p , k) = 



(4.1) 



where we have used the narrow width approximation. To lowest order, the real and imaginary parts of this equation 
become 



r k 



+ U R (u p ,k) = 



-sgn(r fc ) 



nj(^p, fc) 

2 LOr, 



(4.2) 
(4.3) 



and the lowest order solution of eq. (4.2) is given by 
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CJp = ± 



UJk 



U R (uj kl k) 
2uj k 



= ± 



(4.4) 



The solution for the imaginary part is obtained by replacing lu p = iuj k in Hi(u) p ,k) to this order. However, the 



equation for the imaginary part (4.3) does not have a solution because II/(cj, k) is an odd function of w and Hi(u>k, k) > 
0. Therefore, there is no complex pole in the physical sheet. This is a fairly well-known (but seldom noticed) result: 
if the imaginary part of the self-energy on shell is positive there is no complex pole solution in the physical sheet, the 
pole has moved off into the unphysical (second) sheet. 

In the case that Hi(u)k,k) is negative, complex poles appear in the physical sheet, but in such case there are 
two poles with both signs for Tk, i.e., one corresponds to a growing exponential in time and the other a decaying 
exponential, this is the signal of an instability, not of damping. However since we confirm that Hj(ujk,k) is positive 
in the case under consideration, the complex poles are in the unphysical sheet, describing a resonance. 



B. Scalar decay implies fermion damping 

Since we have determined that there are no complex poles in the physical sheet in the s-plane, i.e., the integration 
region, the only singularities are the cut discontinuities along the following segments of the imaginary axis: 



i\Jk 2 + (m + M) 2 , ic 



-i\]k 2 + (to + M) 2 , -ic 



and 



-iy/k 2 + (to - M) 2 , iyjk 2 + {m- Mf 



The contour of integration (Bromwich contour) is deformed to wrap around these cuts. In the narrow width approx- 
imation and consistent with perturbation theory, the discontinuities in the self-energy are perturbatively small; and 
since either the real pole or the resonance is below the two particle cut, the contribution from this cut is always per- 
turbatively small. On the other hand, in the relevant case of m > 2M with a quasiparticle resonance, the contribution 
from the cut 



= -i\]k 2 + (to - M) 2 , iy/k 2 + (m - Mf 



becomes the dominant one. It is convenient to write the product in eq. (2.24) in the simplified form 

Sf(s, k) 



S(s,k) (j-k + M- E pi g(0)) = 



-s 2 ~u>i + n( s , k) ' 



(4.5) 



which defines J\f(s,k), and to change variables to s = —iu> ± + on both sides of the cut, leading to the following 
contribution from the thermal cut to the real time evolution 



N R (uj,k) n/(w,fc) 



, 2 - uj 2 + ii R {oj,k)] 2 + ni(u,ky 
Mi{u,k) [uj 2 -Q 2 + n R (uj,k)} 



n R (uj,k)} 2 + nj(uj,ky 



(4.6) 



where J\f R j(oj, k) are obtained by replacing the el (uj) by their real or imaginary parts and lo^ — ±y/k 2 + {m — M) 2 . 

The term proportional to j\f R (u>, k) features a typical Brcit-Wigner resonance shape near the real part of the complex 
pole at us 2 — uj 2 + H R (ojp, k) = since, for to > 2M, the imaginary part of the self energy at this value of u> (per- 
turbatively close to ±Wfc) is non- vanishing. On the other hand, the term proportional to Mi(lo 1 k) is a representation 
of the principal part in the limit of small Hi(u>,k) and is therefore subleading. The sharply peaked resonances at 
lu = ±|cj p | ~ ±d>fc dominate the spectral density and give the largest contribution to the real time evolution. In the 
limit of a narrow resonance, the u> integral is performed by taking the integration limits to infinity and approximating 
near the resonances at uj — ±|cj n ' 
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n R (uj,k)] 2 + n/(w,A?) 2 2w2( WT | Wp | 



r 2 fe ' 



an fl ( w ,fc) 



Z fc IIj(^p, fc) 
2u> n 



(4.7) 



The wave- function renormalization constant Z k is finite since the self-energy has been rendered finite by an appropriate 
choice of counterterms. Only when the counterterms are chosen to provide a subtraction of the self-energy at the 
position of the resonance uj p will result in Z\. = 1. On the other hand, if the counterterms are chosen to renormalizc 
the theory on-shell at zero temperature, then describes the dressing of the medium and is smaller than one. 

The integration in the variable uj is performed under these approximations (justified for narrow width) leading to 
the real time evolution 



Z k 



"A^(|^ p |,fc)e- l l^l*+A^(-|^ p |,fc) e ^l t ' 



-r k t 



(4.8) 



Gathering the results for the imaginary part of the effective self-energy in eq. (3.9), we find the real-time damping rate 
for the fermionic mean fields to one loop order for arbitrary values of the scalar and fermion masses (with m > 2M), 
temperature and fermion momentum to be given by 



Try 



d 3 q WkUJq —k-q— M 2 
(2tt) 3 2u} k 2Co q uj k+q 



y 2 m 2 T 
32Trfco>i. 



4M< 



In 



(nk+q + n q ) S(uj k +UJq- UJ k+ q) 



1 



(4.9) 



where q* and q% are given by eqs. ( p.lC| ). 

Figures (1-3) display the behavior of F^ for several ranges of the parameters. We have chosen a wide range of 
parameters for the ratios of the scalar to fermion masses (m/M) and the temperature to fermion mass (T/M) to 
illustrate in detail the important differences. The damping rate features a strong peak as a function of the ratio 
k/M. This peak is at very small momentum when the ratio of scalar to fermion mass is not much larger than 2, but 
moves to larger values of the fermion momentum when this ratio is very large. Figure 1 displays this feature in an 
extreme case (m/M = 800) to highlight this behavior. The height of the peak is a monotonically increasing function 
of temperature as expected. 

This is one of the important results of this work: the induced decay of the heavy scalar into fermion pairs results 
in a damping of the amplitude of fermionic excitations. 



C. Resonance wave functions and all-order expression for the damping rate 



Since there is no complex pole solution in the physical sheet, there are no solutions of the effective in-medium Dirac 
equation for m > 2M. However, we can define the spinor wave function of the resonance by considering the solutions 
of the in-medium Dirac equation with only the real part of the self-energy 1 10 ^] at the the value uj — uj p . Since the 
form given by eq. (2.23) for the self-energy is general and not restricted to perturbation theory, our analysis below is 
valid to all orders. 

-iuj p and the real part of the coefficients sv (u> p ) 



Using the expression for the self-energy given by eq. (2.23) for s 



given by eqs. (3^, 3^), we now introduce the following variables 



W P = w P 



■ e (0) -K)] , 

R,k y p;l ' 



M = M\1 



then the resonance wave functions ^(wp) obey the following effective in-medium Dirac equation 



(4.10) 



(4.11) 



Using the property (uj) = e^(—uj) (see eq. (p\4|)), we define the two particle solutions as U ? = ^(^p) f° r 
uj p = \uj p \ and the two antiparticle solutions as V.g = ^ , _r(a; p ) for lo p = — \u> p \. These in-medium solutions satisfy 
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M)U s . s = and (jC + M)V s . % = 



(4.12) 



with |C = 70 1 W p | — 7 • /C. 

These Dirac spinors can be ob tained from the us ual free p a rtic le solutions by the replacement ui k 



and M —> M. given by eq. ( 4.10 ). Combining eqs. ( 2.23 , 3.5 , 3.6 ) and the definition of the width given by (4/7), we 
obtain the following expression for the width of the resonance to all orders in perturbation theory 



4|w p | 



7Tr[ 



■X)S/(|wp|)] 



2|w D 



s=± 



[^.gS/(|wp|)^.g 



(4.13) 



where we have alternatively written the expression for the width in terms of the exact resonance spinors in the medium. 
This result confirms those found in reference K,E8| and leads to the often quoted expression for the width (7||] in 
lowest order. 



V. KINETICS OF FERMION RELAXATION 



To clarify and confirm independently the result of the previous section of induced decay of the scalar leading to 
a quasiparticle width of the fermionic excitations, we now provide an analysis of the relaxation of the distribution 
function for the fermions via a kinetic Boltzmann equation. This analysis will also provide a firm relationship 
between the damping rate and the interaction rate in real time in the relaxation time approximation (linearization 
near equilibrium). 

Let us denote the distribution function for scalars of momentum k at time t by N^(t) and that for fermions 
of momentum k and spin s by Ns (t) . In the kinetic approach, the derivative of this distribution function with 
respect to time is obtained from a Boltzmann equation. Since for a fixed spin component the matrix elements 
for the transition probabilities are rather cumbersome, we define the spin-averaged fermion distribution function as 

Two processes are responsible for the change in the fermion populations in a fermion-scalar plasma: i) <fi — > ip + ip 
(creation) which provides the 'gain' term in the balance equation, and ii) ip + ip — * <p (annihilation) which provides the 
'loss' term. Using the standard approach to obtain kinetic Bolzmann rate equations, we find that the spin-averaged 
rates for creation and annihilation are given by 



d 
df 



d 3 q uik^q — k ■ q— M 2 



d 
It 



gam - (27r) 3 2cJ 9 UJ k UJk+ q 

Nk+M 1 ~ ~ N q -(t)]5(u k +co q - uj k+q ) , 

2 f d 3 q Lb k uj„ — k ■ q— M 2 
— / 7 — 7^-^ x 

loss J (27r)' J 2w g UkU k +q 

[! + N n+ lt)]N % {t)NS)^k +u q - u; k+q ) , 



respectively. The spin-averaged net rate is simply 



IjVgfi) - ^Nr(t) 
dt fcV ; dt feV ' 



gam 



d 
dt 



T,N k (t) 



(5.1) 



loss 



Let us now consider that all of the modes but the fermionic mode with wavevector k in the fermion-scalar plasma are 
in thermal equilibrium, while the population for the fermionic fc-mode has a small deviation from thermal equilibrium, 
i.e., 



1 



+ SNdt) , N ? (t) 



1 



e Pu) k _|_ 2 fc V 1 ' qK ' g/3w 9 _|_ I ' k+q\ I e f3u k+q _ ]_ 

In the linear relaxation approximation (or relaxation time approximation), we find 



(5.2) 



-SNtf) = -T fe SN^t) , 



(5.3) 
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with Tfc being the interaction rate in the relaxation time approximation (linear relaxation) and given by 



y 2 m 2 T ( 4M 2 ', 
t^— I 1 — I In 



(5.4) 



where q\ and q% are given by eqs. (3.1C). Comparing the interaction rate with the damping rate found in the previous 



section (see eq. (4.9)), we provide a real time confirmation of the result 

T k = 2T k . (5.5) 

The kinetic analysis confirms that the damping of the fermionic quasiparticle excitations in the medium is a 
consequence of the induced decay of the heavy scalar. Furthermore this analysis in real time clearly establishes the 
relation between the interaction rate in the relaxation time approximation and the exponential decay of the amplitude 
of the mean field at least to lowest order. Recently a detailed investigation between the damping and interaction rates 
for chiral fermions in gauge theories at finite temperature has been reported in reference p9[ | within the framework of 
imaginary time (Matsubara) finite temperature field theory. Our analysis provides a complementary confirmation of 
this result in real time both for the relaxation of the mean field and that of the spin-averaged distribution function. 

VI. CONCLUSIONS, COMMENTS AND FURTHER QUESTIONS 

In this article we have focused on studying the propagation of fermionic excitations in a fermion-scalar plasma, as 
a complement to the more studied issue of propagation in a gauge plasma. 

Our motivation was to provide a real time analysis of the propagation that could eventually be used in other problems 
such as, for example, neutrino oscillations in medium and in non-equilibrium processes in electroweak baryogencsis. 
The first step of the program is to obtain the effective Dirac equation in medium and in real time. This is achieved by 
relating the problem of linear response to an initial value problem for the mean field that is induced by an external 
Grassmann-valued source term. The resulting Dirac equation for the mean field is fully renormalized and causal and 
allows a direct study of real time phenomena. 

We used this description to study the propagation of fermions in a fermion-scalar plasma to lowest order in the 
Yukawa coupling. We found that when the scalar mass is large enough that its decay into fermion pairs is kinematically 
allowed, this process in the medium leads to a damping of the fermionic excitations and a quasiparticle picture of its 
propagation in the medium. A real time description of this process clearly leads to the identification of the damping 
rate, which we computed to one loop order for arbitrary values of the fermion and scalar masses (provided the scalar 
is heavy enough to decay), temperature and fermion momentum. 

An all-order expression for the damping rate (in the narrow width approximation) is obtained from the exact 



quasiparticle solutions to the in-medium Dirac equation [see eq. ( 4.13 )]. A kinetic approach based on a Boltzmann 
equation for the spin-averaged fermionic distribution function reveals that the interaction rate in the relaxation time 
approximation (linear departures from equilibrium) are simply related to the damping rate of the mean fields at 
least to lowest order in the Yukawa coupling. We emphasize that this relation is established here from the real time 
evolutions both of the mean field and the distribution function. 

Comments. Although the expression for the rate was obtained for arbitrary values of the scalar and fermion masses, 
temperature and fermion momentum, a deeper analysis is required if the theory undergoes a second or very weakly 
first order transition. The reason being that if the fermionic masses are a result of spontaneous symmetry breaking 
in the scalar sector, near a second order (or very weak first order) phase transition both the scalar mass and the 
chiral breaking fermion mass vanish. In this case the kinematic region in momenta for which the energy conserving 
delta functions arc fulfilled shrinks and one must understand if, for soft fermionic momentum, a resummation akin 



to hard thermal loops is required. In particular from the expression of the spectral densities given by eq. ( |3.2| ), it is 
straightforward to see that both p^(ko) and pj? (&o) contribute in the hard thermal loop limit for vanishing scalar 
and fermion mass very similarly to the Landau damping contribution from gauge fields |l7| with a cut discontinuity 
for space like momenta (uj 2 < k 2 ). This particular case, corresponding to the limit T S> to, A/, would have to be 
studied separately and is beyond the scope of this article. A possibility that arises in this situation is that as the phase 
transition is approached the effective temperature dependent scalar mass becomes smaller than twice the temperature 
dependent fermionic mass, the scalar decay channel shuts-off and the contribution to the fermionic damping rate 
vanishes. This of course depends on the self-couplings of the scalar sector and requires studying in detail particular 
models. 

Further questions. In a full fermion, gauge plus scalar theory with light fermions of mass M -C eT (e is the gauge 
coupling), the contribution of the gauge sector to the fermion self-energy requires the hard thermal loop resummation 
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for soft fermion momenta k of order eT, but is perturbative for hard fermion momenta k of order T. For the lightest 
quarks u, d and s but certainly not c, b and t (here we have in mind the electroweak theory and the problem of 
baryogenesis) , the Yukawa couplings y <C e and the contribution from the scalars is perturbatively small. This aspect 
notwithstanding, let us consider the case of soft fermionic momenta. The hard thermal loop resummation leads to 
the dispersion relations for plasmino quasiparticles in the medium, which in lowest order in HTL are stable collective 
excitations. However, with the coupling of scalars, the heavy scalar (with mass of order T) can now decay into a soft 
plasmino and a hard fermion (as the hard plasminos are almost indistinguishable from the zero-temperature fcrmions). 
The results of this article indicate that this process will then lead to a damping rate for the plasminos in lowest order 
(and certainly perturbative) in the Yukawa coupling. This process will compete with the damping of soft fermions 
through the exchange of magnetic photons, and it requires a detailed analysis which is currently under study |p4j| . 
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k/M 

FIG. 2. 32Trr k /(My 2 ) v.s. k/M for m/M = 4 and T/M = 150 (solid line), 100 (dotted line), 50 (dashed line). 
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k/M 

FIG. 3. 32nr k /(My 2 ) v.s. k/M for m/M = 4 and T/M = 1000 (solid line), 800 (dotted line), 600 (dashed line). 
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